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Coalescing binary neutron stars (NS) are expected to be an important source of gravitational 
waves (GW) detectable by laser interferometers. We present here a simple method for determining 
the compactness ratio M/R of NS based on the observed deviation of the GW energy spectrum 
from point-mass behavior at the end of an inspiral event. Our method is based on the properties of 
quasi-equilibrium binary NS sequences and does not require the computation of the full GW signal 
h(t). Combined with the measurement of the NS masses from the GW signal during inspiral, the 
determination of M/R will allow very strong constraints to be placed on the equation of state of 
nuclear matter at high densities. 
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Coalescing compact binaries containing two neutron stars (NS) are among the most important sources of gravita- 
tional waves (GW) for LIGO jj], VIRGO and other laser interferometers. Should the inspiral of such a binary be 
detected, the frequency evolution of the GW signal will immediately yield the system's "chirp mass" M c u = (^"M 2 /". 
where /i and Ai are the reduced and total mass of the binary, respectively. Higher-order post-Newtonian effects on 
the phase evolution of the signal also allow for the determination of the reduced mass fi, and thus the individual 
masses Mi and M2 of the two NS |||. The determination of the NS radii in addition to their masses would yield 
important information about the equation of state (EOS) at nuclear densities, and could even indicate the presence 
of more exotic phases, such as strange quark matter instead of ordinary nuclear matter Q. The GW signal of a 
coalescing binary could yield such information but this is limited in two different ways. During the slow inspiral phase 
at large separations, i.e., low frequencies (/ <C 1kHz), the stars behave like point masses, and finite-size effects are 
not expected to leave any signature in the GW signal J3], ||, During the final hydrodynamic merger, characteristic 
frequencies of GW emission could yield important information about the fluid EOS |7], ^, [|, but these frequencies 
are expected to lie well above the photon shot noise limit of current interferometers (/ > 1kHz). Thus, it is only 
during the last few orbits of the inspiral, just prior to merger, that we can hope to see the imprint of the NS radii on 
a measurable GW signal (with / < 1kHz). 

Several groups have studied this terminal phase of inspiral by constructing quasi-equilibrium sequences of close NS 
binaries in the conformal flatness approximation of general relativity (GR) (l^, [ll], 12], In this approximation, 
it is assumed that the binary system evolves along a sequence of appropriately constructed equilibrium states with 
decreasing binary separation as energy is radiated away in GW. From the binary equilibrium energy curve E cq (r), 
which gives the total system energy as a function of binary separation r, and the GW luminosity Lew, one can 
derive the radial infall rate as v r = LQ^(dE cq / dr)~ l . With v r — dr/dt known, this provides the time evolution along 
the equilibrium sequence and the GW signal h(t). This approach should remain accurate as long as the radial infall 
timescale r/v r is longer than the dynamical timescale of the system, i.e., until the point where dynamical instability 
sets in, and the two stars plunge inward rapidly and merge. 

Unfortunately, calculating the correct GW luminosity for a given matter configuration in GR is an extremely difficult 
task. Different approaches have required either time integration of the full non- linear equations of GR H], or the 
solution of a complicated wave equation for terms representing the spherical harmonic expansion of the GW metric 
perturbation ]T^| . The great complexity of these approaches is in stark contrast with the simplicity of the quasi- 
equilibrium approximation. However, we point out here that the GW energy spectrum dEcp^/df can be calculated 
directly and very simply from the equilibrium energy curve, independent of any knowledge about the GW luminosity. 
Indeed, by definition of the quasi-equilibrium approximation, the energy decrease —dE cq between two neighboring 
binary configurations along the sequence is equal to the energy dE'cw radiated away as the wave frequency sweeps 
up by df, where the GW frequency is twice the orbital frequency, / — 2/ or b- Thus, one should simply compute the 
total energy E cq as a function of frequency / along the equilibrium sequence, and the GW energy spectrum is then 
given by the derivative dEsw/df — —dE cq /df. As a trivial example, consider the inspiral of two point masses in the 
Newtonian limit, where we have E cq oc — and / cx r -3 / 2 . It follows that E cq oc — J 2 / 3 and thus dEow/df cx f^ 1 ^ 3 , 
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a well-known result. In addition to the assumptions underlying the quasi-equilibrium approximation, the validity of 
this simple approach relies on the additional assumption that the GW emission during the later merger phase has no 
effect on the energy spectrum at lower frequencies. Indeed, this has been demonstrated in numerical hydrodynamic 
calculations of binary mergers || , which show a clear separation between the inspiral and merger components of 
the emission in frequency space. 

We have investigated the properties of the GW emission during the final phase of binary NS inspiral using new, highly 



accurate equilibrium sequences calculated with the Meudon code [16 . The formalism is based on the conformally flat 
approximation, which is expected to yield accurate equilibrium matter configurations for this phase [[l0| [ll], Jl2, 13 . 
The resulting five non- linear, coupled elliptic equations are solved using a multi-domain spectral method [17, 18 1. 
This approach has already been used successfully in various astrophysical applications jl3L |l9], ^(j. Typically, the 
computed fields satisfy the constraints of full GR to within ~ 1% j2lj. The code has been improved recently, especially 
with regard to the treatment of the external compactified domain, and the numerical error in computing equilibrium 
configurations, measured in terms of how well the virial theorem is satisfied, is of order one part in 10 5 . The equilibrium 
sequences presented here are the natural extension of the work already published in Jl3|, |2^] . They will be discussed 
in more detail in ^3|. Here we only show the variation of the ADM mass of the system (total binary mass-energy 
M. = c~ 2 E oq ) and the GW frequency (twice the orbital frequency), which are sufficient for our purposes to determine 
the GW energy spectrum. 

The equilibrium configurations have been calculated for irrotational binaries, i.e., assuming that the fluid has zero 
vorticity in the inertial frame. Indeed, the NS should be spinning slowly at large separations and the viscosity of NS 
matter is too small for tidal spin-up to become significant on the coalescence timescale fi |24|] . Based on the current 
set of well-measured NS masses in relativistic binary radio pulsars, it is expected that all NS in coalescing binaries will 
have masses M ~ 1.35M Q |25j] . Hence, for simplicity, we consider only equal-mass binaries where Mi = M 2 = 1.35M Q . 
Also for simplicity, we model the NS EOS with a simple polytropic form P = Kp T , where P is the pressure and p the 
rest-mass density. The constant K represents the overall compressibility of the matter and largely sets the value of 
the stellar radius for a given mass, while the adiabatic exponent T measures the stiffness of the EOS and affects the 
degree of central concentration of the NS interior. Based on our experience with hydrodynamic calculations we 
expect that the GW energy spectrum just prior to merger is determined primarily by the NS radius through K, with 
relatively little sensitivity to T. For this reason, in this initial study, we allow K, and therefore also the stellar radius 
R, to vary for different NS models, but we set T = 2 for all models, as this value fits well most published NS EOS 
(see, e.g., [^6| and references therein). Specifically, we consider NS models with compactness ratios M/R =0.12, 0.14, 
0.16, and 0.18 (setting G = c = 1), where M is the ADM (gravitational) mass measured by an observer at infinity 
for a single isolated NS, and R is the circumferential radius of the NS. For M = 1.35 Mq, the corresponding radii 
are R = 16.6, 14.2, 12.4, and 11.1 km, respectively, spanning the range of values for NS radii calculated from various 
physical EOS. Note, however, that our results are unchanged under the rescaling given by R' = kR, M' = kM, 
f = f /k, for any constant k. 

For each NS model, about 12 binary equilibrium configurations are computed with decreasing separations, until a 
cusp develops on the NS surface. Equilibrium configurations for smaller separations do not exist. To each sequence 
we fit a curve of the form 

M(f) = 2.7M G -k N f/ 3 + k 1 f + k 2 f 2 (1) 

to represent the variation of total mass-energy as a function of GW frequency. The first term gives the total gravita- 
tional mass of the system at infinite separation, while the second term represents the Newtonian point- mass behavior, 
with k N = 2" 4 / 3 7r 2 / 3 G 2/3 M 5 / 3 = 4.056 x 10~ 4 M Q Hz~ 2/3 . The term cx / was introduced heuristically to represent 
the lowest-order post-Newtonian and finite-size corrections to the point-mass behavior at intermediate frequencies. 
The term oc f 2 represents the tidal interaction energy, which causes the equilibrium energy curve to flatten at high 
frequencies. Our best (least-squares) fit values of k\ and k 2 for each sequence are listed in Table 1 and the results 
are illustrated in Fig. [l]. The asterisks show the data points along each sequence, with a typical error between the 
data points and the fit of 5M ~ 1O~ 4 M . We find in all cases that k 2 is positive, as we would expect: tidal de- 
formations and relativistic gravitational effects increase the equilibrium energy p7| , p8| . We note that none of the 
equilibrium curves shows evidence of an energy minimum, which would have implied the onset of dynamical instability 

Computing the GW energy spectrum for each model now only requires differentiating the fitted curves with respect 
to frequency. The results are shown in Fig. |^. In each case, we see a characteristic frequency range where the spectrum 
plunges rapidly below the extrapolation of the low- frequency result. This corresponds to the flattening of the energy 
curves in response to the growing tidal interaction and PN effects. Also shown is the energy spectrum of a 3PN, 
irrotational, point-mass binary, computed according to the results found in Ref. [BOl, which closely tracks our most 
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compact model, indicating that the differences we see in the energy spectra result from finite-size effects associated 
with the NS radius. To quantify the deviations from the Newtonian case, we define a set of break frequencies, at which 
the energy spectrum has dropped by some factor below the point-mass result. The values we find for / 10 , /2s, and 
/ 50 , where dE GW /df has dropped by 10%, 25%, and 50%, respectively, are listed in the last three columns of Table 1. 
We see that all these characteristic frequencies lie within the frequency range accessible by LIGO-type detectors, 
with perhaps the exception of /50 for the more compact sequences. Note that the calculation of /50 values requires 
extrapolating the equilibrium energy curves slightly beyond the last equilibrium model (where a cusp develops), and 
may therefore be less reliable. However, we find that /50 has a particularly steep, quasi-linear dependence on the NS 
compactness, given by /50 — [10 4 (M/i?) — 460] Hz within the range of NS radii we considered. For comparison, ^25 
values can be determined safely within the quasi-equilibrium approximation, and the sensitivity on compactness is 
only slightly reduced, with / 25 ~ [5000(M/i?) - 85] Hz. 

The best definition of the break frequency will be a tradeoff between higher signal-to-noise ratio at lower frequencies 
jnj, and greater ease of discrimination between different EOS at the higher frequencies. In addition, the quasi- 
equilibrium approximation is expected to be most accurate at lower frequencies, where the inspiral rate is slower. 
However, we doubt that this could become a major issue: if we adopt, for simplicity, the point-mass formula for 
the GW luminosity, and compute the corresponding radial infall velocity along our equilibrium sequences, we find 
that v r never exceeds 5% of the orbital velocity, even at the point where we define /50 (the corresponding fraction 
at the point where we define /2s is about 2%) j32|. Ultimately, the optimum choice should be determined by data 
analysts, taking into account the accuracy with which the break frequencies can be extracted using matched filtering 
or other techniques |3J| . Preliminary studies of this problem have already been performed for both broad-band and 
narrow-band interferometer configurations |34|, Hi- Defining a precise break frequency may not even be necessary. 
Instead, the GW inspiral templates could be terminated at high frequency in a manner that reproduces the energy 
spectrum given by a simple analytic form, such as our Eq. [l]. The free parameters k\ and fc 2 could then be measured 
experimentally and compared directly to the predictions of various realistic nuclear EOS used in computing binary 
equilibrium sequences. In future work, we plan to compute such sequences, and the corresponding energy curves, for 
a wide variety of published, realistic NS EOS. 

We are grateful to Kip Thorne for originally pointing out the importance of the break frequencies in GW spectra. 
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TABLE I: Properties of the quasi-equilibrium sequences. Here M/R (with G = c = 1) is the compactness of an isolated NS 
seen by an observer at infinity, R is the circumferential radius for an ADM mass of M = 1.35 Mq, / c is the GW frequency 
at the final point of each sequence (cusp), k\ and k<2 are the best fit parameters in Eq. jj], and /io, /2s, and /50 are the break 
frequencies at which the GW energy spectrum has dropped, respectively, by 10%, 25%, and 50% below the point-mass energy 
spectrum. 



M/R 


R (km) 


fc (Hz) 


fcl &2 


fio (Hz) 


/25 (Hz) 


/50 (Hz) 


0.12 


16.6 


641 


-4.939E-6 1.290E-8 
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518 


764 


0.14 


14.2 


807 


-3.363E-6 9.244E-9 
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0.16 


12.4 


1002 


-1.806E-6 6.490E-9 
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720 


1137 


0.18 


11.1 
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-5.834E-7 4.835E-9 


431 


810 


1331 




2.68 
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FIG. 1: ADM mass (total mass-energy) of a binary NS system as a function of GW frequency (twice the orbital frequency), 
computed along each of our 4 irrotational equilibrium sequences. From bottom to top, the sequences correspond to NS with 
compactness M/R — 0.12, 0.14, 0.16, and 0.18. All models assume a polytropic EOS with T = 2 and a NS mass of 1.35 Mq for 
both components. The asterisks indicate the individual equilibrium configurations calculated along each sequence, while the 
lines show our best fit using Eq. [I] and the values of Table 1. 
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FIG. 2: Energy spectrum dEaw/df of GW emission emitted along each of the 4 sequences of Fig. |l| in the quasi-equilibrium 
approximation. Also shown is an irrotational 3PN point mass binary from Ref. [ |30| , which closely tracks our most compact 
model. Asterisks indicate the terminal point along each sequence, where a cusp develops. The slanted straight lines show, 
from top to bottom, the point-mass Newtonian energy spectrum (oc / — 1 ^ 3 ) multiplied by 1.0, 0.9, 0.75, and 0.5. The last three 
values are used to define characteristic break frequencies /10, /2s, and /go, where the energy spectrum has dropped by the 
corresponding fraction. The units on the right and top axes show the corresponding dimensionless quantities, with the mass 
dependence scaled away. 



